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Abstract 

Two new sum rules were recently discovered by Le Yaouanc et al. by applying the operator 
product expansion to the nonforward matrix element of a time-ordered product of b —* c currents in 
the heavy-quark limit of QCD. They lead to the constraints a 2 > 5p 2 /4 and a 2 > 3( / o 2 ) 2 /5 + 4 / o 2 /5 
on the curvature of the B — > Isgur-Wise function, both of which imply the absolute lower 
bound a 2 > 15/16 when combined with the Uraltsev bound p 2 > 3/4 on the slope. This paper 
calculates order a s corrections to these bounds, increasing the accuracy of the resultant constraints 
on the physical form factors. The latter may have implications for the determination of \V c b\ from 
exclusive semileptonic B meson decays. 
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I. INTRODUCTION 

Heavy quark effective theory (HQET) 1] provides a model-independent method of ex- 
tracting the CKM matrix element \V c b\ from exclusive semileptonic B meson decays. The 
B — > D(*Hv differential decay rates are given by 
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D*lu) = g|l H m| r. 3 (l - r.) 2 (w + if 



dw y 487r 3 

4u> 1 — 2wr* + r 2 
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^(B^DW) = ^^r 3 (l + r)> 2 -l) 3 /^ D H 2 , (1) 

where rv*) = m D (*)/mB and w = t> • v ' is the product of the B and four-velocities. 
Heavy quark symmetry [2j relates i? — > D^W form factors to the corresponding Isgur-Wise 
function, with the result J r n*(w) = J-n(w) = £(w) in the heavy-quark limit of QCD. Since 

nnnn 

£(w) is absolutely normalized to unity at zero recoil (i.e., w = 1) |2J,|3|,II|5(, experimental 
data determine \V c b\ without recourse to model-specific assumptions. 

This procedure has several sources of uncertainty. First, the identity ^(^(l) = 1 receives 
both perturbative corrections and corrections suppressed by the heavy b and c quark masses. 
The former are known to order az with unknown higher-order corrections likely less than 
1%, but the latter depend on four subleading Isgur-Wise functions and have been estimated 
only with phenomenological models and quenched lattice QCD. 

Another source of error is the extrapolation of measured form factors to zero recoil, where 
the rates vanish. Linear fits of J : D {*){w) underestimate the zero-recoil value by about 3%, 
an effect mostly due to the curvature [2J. Using non-linear shapes for F D {*)(w) reduces this 
error, and therefore constraints on second and higher derivatives at zero-recoil are welcome. 
Dispersive constraints flQ relate second and sondes hi g her derivatives to the first and 
are commonly used. 

HQET sum rules provide a complementary way of constraining the J^w (w) shapes. New 
sum rules for the curvature and higher derivatives of the Isgur-Wise function were derived 
in Refs. 10, 11,^^. Equating the result of inserting a complete set of intermediate states 
in the nonforward matrix element of a time-ordered product of HQET currents with the 
operator product expansion (OPE) gives a generic sum rule depending on the products of 
the velocities of the initial, final, and intermediate states. These are denoted respectively 
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by Vi, Vf, and v' (the intermediate states all have the same velocity v' in the infinite-mass 
limit), and the products are denoted by 



Wif = Vi-Vf , Wi = Vi-V' , Wf=VfV, 



or generically w x . These parameters are constrained to lie within the range 



(2) 



Wi,Wf,Wif >1 , WiWf - \J (wf - l){wj - 1) < w if < WiW f + yj (wf - l)(wj - 1) , (3) 

and differentiating the generic sum rule with respect to them at w x — 1 (read: W{ = Wf = 
Wif = 1) produces a class of sum rules for derivatives of the Isgur-Wise function at zero 



recoil. The sum rules of Refs. 



10 



11 



12 1 were derived at tree level in the heavy-quark limit. 



The present paper includes the order a s corrections to the new sum rules and uses them 
to derive bounds on the curvatures of T D (,){w) including ot s and Aqcd/"t. c ,6 corrections. 
Including these corrections to the heavy-quark limit is important for meaningful comparison 
with data and dispersive constraints. 



II. DERIVATION OF THE GENERIC SUM RULE 

The derivation of the generic sum rule follows a well-known formalism 

Q,H HQ- 
it begins with the consideration of the time-ordered product of two arbitrary heavy-heavy 
currents 

Tfi(e) = ^fd 4 x (B f (p f )\T{J f (0), J t (x)} (4) 



as a complex function of e — Em — Ei — q° at fixed q, where Em = \J m2 M + \Pi + Q I 2 is the 
minimum possible energy of the charmed hadronic state that Jj can create at fixed q. The 
currents have the form 

J f (x) = b(x)T f c(x) , Ji(y) = c(y)T t b(y) (5) 

for any Dirac matrices Tij. Only the choices Tij — and Tij = fafjs are explored 
here, but in principle others lead to different sum rules. The B states are ground state 
B or B* mesons and have the standard relativistic normalization. As in the derivation 



of the Uraltsev sum rule 



In l 



the initial and final states do not necessarily have the same 
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velocity. Considering the nonforward matrix element of the time-ordered product is a crucial 
generalization in deriving the new sum rules [lp| . 

From Eq. (J3J) one proceeds by splitting up the two time-orderings and inserting complete 
sets of intermediate charm states. The result is 

1 Vf^WiT ff+3 ) ( B f\M0)\X 5bb )(X 5bb \J f (0)\B t ) 



where the sums include phase space factors such as d 3 p/ (2n) 3 2Ex- Again, T/» has been 
written as above to call attention to the full generality possible for deriving sum rules by 
this method, but here both Bi and Bf will be taken to be B mesons to avoid the considerable 
complication of the B* polarization. In addition, HQET states and currents will be used 
henceforth since the goal is sum rules for the derivatives of the Isgur-Wise function. Deriving 
the bounds in the effective theory also makes the calculation of perturbative corrections much 
easier. 

The function Tfi(e) has two cuts along the real shown in Fig. U The important one 

here, running from — oo to the origin, comes from the first time-ordering and corresponds 
to intermediate states with a c quark or a c quark, a b quark, and a b quark. The cut 
associated with the other time-ordering begins near 2m c and corresponds to states with two 
b quarks and a c quark. Since Tfi(e) is perturbatively calculable only when smeared over a 
large enough range of e , it is multiplied by a weight function Wa and integrated around 
the contour shown in the figure. The scale A gives the extent of the smearing and therefore 
should be well above Aqcd- The contour chosen eliminates all but the intermediate states 
with heavy quark content c by avoiding the second cut and pinching the first at e = —Im^. 
Crossing the contour assumes local duality at the scale m&, but if A < m b the weight function 
will be quite small here. This will be clear with the specific weight function used below. 
Assuming it is analytic in the shaded region of Fig. ^ the result is 

^ / «fe W A (e) T li{ e) = E W A (E M -E xJ ^MWq^mmi , (7) 

where the delta function has been used to perform the phase-space integral and the HQET 
state normalization convention has been used to eliminate mass factors in the denominator. 
The intermediate X c states carry four-momentum p Xc = m Xc v' = Pi + q. 



A 



C 



FIG. 1: The cuts of Tfi in the complex e plane. The depicted contour picks up only contributions 
from the left-hand cut, which corresponds to physical states with a charm quark. The states given 
by the right-hand cut do not contribute here. 
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Choice of the weight function is governed by well-known concerns 
one uses W&(e) = 6 (A + e), which is the n — > oo limit of the set of functions 

A 2n 



Q. 



w£\e) 



In practice 



(8) 



£-2n _j_ /\2n 

for e < 0. But since the weight function must be analytic within the contour, the use of these 
is strictly correct only for small n > 1. In this case the poles at e — 2 \/—l A are a distance of 
order A away from the cut, and the contour can be deformed away without getting too close 
to the cut and relying on local duality at a scale below m^. This is not true of the n — > oo 
limit, in which the poles approach the cut and the contour must pinch the cut at the scale 
A instead of m&. This is a problem because the contribution at A is weighted much more 
heavily than that at m/,, and thus the results will depend more strongly on the assumption 
of local duality. However, an explicit calculation shows that the results here do not depend 
on n, just as the authors of Ref. 3| found in their derivation of corrections to the Bjorken 
sum rule. This is not true in other cases, such as the Voloshin sum rule [if]]. The weight 
function in what follows can therefore be considered a simple step function excluding states 
with excitation energies greater than A. Although increasing A includes more states and 
weakens the bounds, the cutoff energy must be chosen large enough to make perturbative 
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FIG. 2: Diagrams contributing to the order a s corrections to the sum rules. The squares indicate 
insertions of the currents Jj and Jf, respectively. The current Jj inserts momentum q, while the 
current Jf carries away momentum q' sufficient to leave the final 6-quark with velocity Vf. The 
velocity-labeled quark fields are those of the heavy quark effective theory. 

QCD appropriate. Choosing A > 2 GeV should therefore be sufficient. 

The sum rule is derived by performing an operator product expansion on the time-ordered 
product of currents on the left-hand side of Eq. (J7J) while writing out the right-hand side 
explicitly in terms of excited-state Isgur-Wise functions. The leading-order OPE relevant for 
B matrix elements consists of a single dimension-three operator, h Vf Tf{l +rf')Tib Vi . Higher 
dimension operators will be neglected here, as they give corrections suppressed by powers 
of Aq CD /A or AQ CD /m c fc . The order a s corrections to the Wilson coefficient of the leading 
operator are given by a matching calculation involving the diagrams in Figs. El and EJ The 
generic sum rule resulting from this is 



where the function F contains the one-loop corrections. In principle, F could be defined 
to include perturbative corrections of all orders. The form of the corrected sum rule would 
be the same since HQET vertices are spin independent. The right-hand side is written out 
explicitly in the next section. 

Working in the rest frame of the intermediate hadrons (i.e., v'° = 1) and using the MS 
scheme with dimensional regularization and a finite gluon mass m, the contributions to a s F 
of the graphs in Figs. Et a )H2fd) are, respectively, 



y(w if ) [1 + a s F(w t , w f , w if )} Tr [(1 + j f ) T f {l + j ') ^(1 + ^) 
= Y.^{E M -E Xc ){B{v f )\J f \X c {v')){X c {v')\Jmv l )) , 



(9) 




(10) 
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FIG. 3: One-loop renormalization of the leading operator in the operator product expansion of T${. 
The blob indicates an insertion of this operator, b Vf Tf(l + if)')Tib Vi . The external lines are bottom 
quarks in the heavy quark effective theory. 



2a s ln (^ + \M _1 ) , 4A 2 /-i 21n:r-ln[l + 2a;(l-z)K-l)] . 

S - Wi { —± — . In— =- + / dx [ —- y —. ^ ^ ) , (11) 



37r Jw 1 - 1 V 2 Jq l + 2x(l-x)(wi - 1) 



2as -w f \ - l w / + V w /- 1 J ln l^! + f ^lnx-ln^xq-xX^-l)] | ( |2) 



3tt 7 J /i 2 l + 2x(l-x)(w / -l) 



4a s 

-Wi 



_. t f 1 dx dydz5(x + y + z- 1) ^ ^ m /A) (13) 
3tt J io ayfz 2 -am?/A 2 

where a = 1 + 2xy(wif — 1) + 2x^ (wj — 1) + 2yz{wf — 1), and a s is evaluated at subtraction 
point /i. The contribution of Fig. I2jd) cannot easily be simplified further, but this is no 
limitation since the sum rules require only the first few terms of F in an expansion about 
w x = 1. The graph in Fig. El contributes with a minus sign to the matching calculation for 
the Wilson coefficient, since it gives the renormalization of the leading operator in the OPE, 
and so its contribution to a s F is 

^_|^ lnK+ ^ ln ^. (14) 

This infrared divergence cancels that of the graph in Fig. Efd), leaving a s F independent of 
the regulating gluon mass. 

From the results above, it is not hard to show that F(l,w,w) = F(w, l,w) = 0. This 
important characteristic of the perturbative corrections is true to all orders in a s , as can 
easily be seen. In the limits Vi = v' and Vf = v', one of the currents in the time-ordered 
product is a conserved current associated with heavy quark flavor symmetry and its matrix 
elements receive no perturbative corrections. Because HQET loop graphs do not change 
the matrix structures of inserted operators, perturbative corrections to matrix elements of 
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the other current cancel those of the leading operator in the OPE. Therefore, the function 
analogous to F including perturbative corrections of all orders will still vanish in these limits. 

The sum rules derived here are primarily of interest near zero recoil, making it convenient 
to expand F about w x — 1 with the definitions 

a a F(wi,w f ,Wif) = b ifl (wi f - 1) + bn(wi - 1) + b fl {w f - 1) + ~bif 2 (w if - l) 2 

X i 

+ -b i2 ( Wi - if + -b f2 {w f - l) 2 + b i2f ( Wi - 1)K - 1) 

+ b i2if (wi - l)(w if - 1) + b f2if (w f - l){w if - 1) + • • ■ . (15) 

There is no zeroth-order term since F(l, 1, 1) = 0. This follows from the identities 
F(l, w, w) = F(w, 1, w) = 0, which also imply 

bfi + b if i = ba + b ifl = , (16) 
6/2 + 2b f2if + b if2 = b i2 + 2b i2if + b if2 = . (17) 

These relations between derivatives of the perturbative corrections can be checked at order 
a s with the explicit values 

4A 2 \ 



Hfi - -Vfi 



bfi = -bu 



b®if - b f2if 
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f> 2 



b l2 = 6 /2 = i£ [ i - ln^-^ , 
hf2 = -^l-^ + ln^- 



5 fi 2 



The derivatives above are all specific to the rest frame of intermediate hadrons. This is 
the frame used henceforth. In other frames (e.g., v® = 1) the weight function depends on 
the w x parameters, and the sum rules are more complicated but not qualitatively different. 

III. VECTOR AND AXIAL VECTOR SUM RULES 

When specific matrices Tij are chosen, the generic sum rule in Eq. Q can be written out 
explicitly in terms of excited-state Isgur-Wise functions using Falk's description of HQET 



S 



states of arbitrary spin [l^. The choice T it f =tf i f yields 



£{w if ) [1 + a s F(wi, Wf, w if )] (1 + Wi + Wf + Wif) 



+ l){w f + 1)53 7^I^ S t{ w ii w fi w if)Ys T f+^ 

(=0 «T n 



+ £ fcK «*/)£ ^(^Ori^K)W A (m M -mf ^ 2 ) , (19) 

where the weight function now bounds excitation mass because v'° = 1. The functions 
r i±i/2( u ') are ^ — *■ ^£±i/2 Isgur-Wise functions, where £±1/2 is the spin of the light degrees 
of freedom, (— l) e+1 is the parity of the state, and the label n counts "radial excitations." 
This name is inspired by the nonrelativistic constituent quark model, but these states can 
be anything carrying the other quantum numbers, including continuum contributions, for 
which n would be a continuous parameter and the sums integrals. In that case, the functions 
T e±i/l( w ) w °uld not be Isgur-Wise functions but other B decay form factors. This possibility 
will be downplayed here, with the assumption that such contributions are small in the bounds 
derived here. Experimental input on B — > DnW, for example, is needed to evaluate this 
assumption. 

The quark model also offers an interpretation of £ as the orbital angular momentum 
between the light antiquark and the heavy quark. The relation of this notation to that of 
Isgur and Wise jstj for the lower values of £ is given by 

4% {n) (w)=e n \w) , t[^\ w ) = 2t[%{w) , T$ n \w) = V3r$( w ). (20) 

The function Si takes into account the polarization of an intermediate state with integral 
spin £ and is defined as 

S t = v iv< ■ ■ v Wl v f ^- ■ ■ v flH £ £ (A)"^ £ S" W ' ( 21 ) 

A 

where e^ 1 '"^ is the polarization tensor of the intermediate state. The sum runs over the 
2£ + 1 polarizations. This quantity was reduced in Ref. [icj to the relatively simple form 

1/2 

S e { Wi ,w f ,w if ) = E<^K 2 - " l) k {w iWf - w if Y~ 2k , (22) 

k=0 



with the coefficient 



,(£!) 2 (2£-2k)\ 



Ce ' k ~ { 1} (2£)\ k\(£ - k)\(£ - 2k)\ ■ (23) 
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Using this formula it is easy to show that Eq. (JT5j) reduces to 

2(1 + w)£(w)[l + a.F(l, w,w)\ = 2(1 + w)£(w) (24) 

in the limit Vi = v', confirming that F(l,w,w) = to all orders. The limit Vf = v' gives 
F(w,l,w) = 0. 

The axial sum rule (i.e., Eq. © with T^j = ^^75) can be written out explicitly in the 



same way 



III: 



£(wif) [1 + a s F(wi, Wf, w if )} (-1 + Wi + w f - w if 



+ Se+i(wi, w f , w if )J2 ^i72( w <) r i+i72( w /)^A(m M -m2^) . (25) 



As in the vector sum rule, the masses of the intermediate states are denoted by Tnf±$ 2 . The 
limits Vi = v' and Vf — v' are trivial for the axial sum rule. 

The doublets with spin of the light degrees of freedom si = I + 1/2 and si = I — 1/2 
contain states with angular momentum £, £ + 1 and £ — 1, £, respectively. But only one 
member of each doublet contributes to the sum rules in Eqs. (|T§|) and (J23j) because of the 
choice of currents. This explains the appearance of only one polarization function for each 
doublet in the vector and axial vector sum rules. 

The zero-recoil normalization of the B — > D^*> Isgur-Wise function allows one to write 

2 

£H = 1 - p\w - 1) + y (w -If-... . (26) 

The axial and vector sum rules (i.e., Eqs. ()19j) and (j25|l ) give expressions for p 2 , a 2 , and 
higher derivatives of when differentiated with respect to the parameters w x at w x = 1. 
Different combinations of derivatives yield different relations. In the v'° = 1 frame, the 
sum rules are invariant under the interchange of Wi and Wf, and it is therefore sufficiently 
general to consider only derivatives with respect to Wif and w = Wi = Wf. Because of this 
simplification, this paper only uses derivatives of the vector and axial sum rules of the sort 

QP+1 



dwfjdw q 



(27) 

Wif=W=l 



Derivatives of the vector sum rule with p + q = 2 give expressions for a 2 , while the extra 
factors of (w x — 1) in the axial rule require p + q = 3 for curvature relations. 
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As an illustration of the method, one can easily derive the Bjorken |2fl |21| and Uralt- 
sev ^| sum rules with order a s corrections. For this it is only necessary to consider p+q = 1 
in the vector rule and p + q = 2 in the axial rule. Taking the vector sum rule first, the equa- 
tion given by the p — 0, q — 1 derivative is trivial, but p — 1, q — gives the Bjorken sum 
rule with one- loop corrections 

^ 4 + £ (I - ta f ) + % K«] 2 + ? h»f • < 28 > 

This equation has been written in the familiar notation of Isgur and Wise using Eq. (f2*U|) . 
The upper limit A on the sums stands for a factor of the weight function {trm ~ m x c ) > 
which serves to cut off the sums. Without it the results are divergent, as can be seen by 
attempting to take the A — ► oo limit in the order a s corrections. Note that he subtraction- 
point dependence is the same on both sides of the equation, since Isgur- Wise functions are 
independent of /i at zero recoil while their slopes depend on it logarithmically 22|. This 



equation should be evaluated near // = A to avoid large logarithms in the perturbative 
expansion. 

The lower bound resulting from ignoring the sums in Eq. (12811 is similar to one derived in 

n n 

Ref. 23] but somewhat weaker. As discussed in Ref. |16], this is the result of using different 
weight functions. That of Ref. j^J is effectively given by the phase space of b decay and 
falls off faster with e, thus reducing the contribution of the intermediate states to the sum 
rule and strengthening the resultant lower bound. A similar effect could be achieved here 
by using a smaller value for A, but as discussed above this makes the use of perturbative 
QCD less reliable. 

The p = 0, q = 2 derivative of the axial equation also gives the Bjorken sum rule. The 
p = 2, q = and p = 1, q = 1 derivatives give the same result, which, when combined with 
the Bjorken rule, gives the traditional form of the Uraltsev sum rule: 

E [t$(1)] 2 - E H/1(l)f = 7 - hn - kba + 6/x) = \ , (29) 

n n ^ Z <i 

where Eq. (|2U|) has again been used. This equation receives no unsuppressed perturbative 
corrections. (There are in fact perturbative corrections suppressed by Aq CD /A but 
such corrections are being neglected here.) In this particular derivation of the Uraltsev 
rule, this is the result of the relation in Eq. (fTT)|) between the first derivatives of F. But 
another derivation from different derivatives of the axial and vector sum rules gives a s 

11 



corrections proportional to F(l, 1, 1) = 0. It appears that the Uraltsev rule is always 
protected from perturbative corrections by the general identities F(l, w, w) = F(w, 1, w) = 
0. This convergent sum rule indicates that Ty\{l) and 73/2(1) become equal as n — > oo. 
Combined with Eq. (J28|) . the Uraltsev rule improves the Bjorken bound significantly: 



pV) = ! + M/^)+3£K1(i: 



3 4a s 
> - + — - 

4 9tt 



3 p? , 



(30) 



Because the Uraltsev rule is not corrected, the corrections to this improved bound are just 
those of the original Bjorken bound. In particular, they are not substantially increased, as 
one might expect from the drastic improvement to the bound. 

Constraints on the curvature of the Isgur-Wise function are obtained from higher deriva- 
tives of the same equations. The three second derivatives of the vector equation and the 
four third derivatives of the axial can be reduced to five linearly independent relations, as 
demonstrated in Ref. Q). Complete with the one- loop corrections derived here, they are 



P 



a' 



a' 



a' 



a 



4^ (l)(n) m (l)(n)' m 3^ (l)(n) nw (l)(n)'^ 2 
7L T 3/2 l-U T 3/2 + r L r l 



E 



T, 



3/2 1 



Dr. 



(i)W 



(2)(n) 
5/2 1 

5 A 



3/2 
2 



l)+b ifl p 2 
+ 2b ifl p 2 - b if2 



2 (l) T l/2 
1 



i) + nhn , 



i/2 



(b 



i2if 



5 

bf 2 if 



^ n 

3 Ar mvnV^xi 2 4, o 4. 6 



5 

2b ifl p 2 - 6 </2 - -b if i ■ 



(31) 
(32) 
(33) 
(34) 



+ -6i/ip - -6</i - - + 6/m/) 



- -6</ 2 - ^ ( 6 a + 26 i2/ + b f2 ) 



(35) 



The last two equations give the bounds of Ref. [12], complete with order ot s corrections. Only 
a couple of orbital excitations occur and in positive-definite form, allowing the derivation of 
nontrivial lower bounds. Using the values of Eqs. (|18|) gives 

193a, 



32a s ln 4A 2> 
457r y? 
20a, 



6757T 457r 



13a,, 4A 2 
In — — 



4a Sl 4A r 
m — — 

9tt 



148a s Aa s , 4A 2 
— — - + — — In— 

6757T 457T 



(36) 
(37) 



where the identity r}?^ ' 1 ( 1 



/i" y b/D7T 407T /i z 

x / 2 v- L ; — ~P 2 nas been used. As demonstrated below in the derivation 
of physical bounds, the subtraction-point dependence is the same on both sides of these 
inequalities. 
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IV. PHYSICAL BOUNDS 



When combined with a s and Aqcd / ^c,b corrections from matching HQET onto the full 
theory, the curvature bounds derived above imply bounds on the zero-recoil derivatives of 
the functions J-' D ^){w). It is convenient to expand these functions about the zero-recoil point 
according to 



^w(l) 



(38) 



A simple matching calculation, taken from Ref. |22j, yields the relations between the Isgur- 
Wise derivatives and those of the physical shape functions: 

4a., , ml a s ( 20\ A 



2 

Pd(*) 



2 

17 D(*) 



7T 



a 2 (p) + 2p 2 (p) 



4a s m\ 

In— £ + 

9tt p? 



27/ 2m, 

C 

dm 27 



°D(*) ' 



7T 



+ 



4a s m 2 
In— £ 

157T /i 2 



, a * f A («.) 16 ^l , A a(V«0 



7T 



(39) 



The perturbative corrections are model independent. The parameters S^"J } and A^/j are 
given by 

2(1 - z)(ll + 2^ + llr 2 ) + 24(2 -z + z 2 )z \a.z 



>D* 



'D 



27(1 -^) 3 

2(1 - z)(23 - 34z + 23z 2 ) + 12(3 - 3z + 2z 2 )zlnz 
27(1 - z) 3 
8(47 + 17z + 252z 2 + 17z 3 + 47z 4 ) 



0.24 , 
1.20 



(«.) 



675(1 - z) 4 
4(5 + 125z - 55z 2 + 95z 3 - 18z 4 )z\nz 
135(1 - zf 
4(47 - 258z + 302z 2 - 258z 3 + 47z 4 ) 



-1.16 



225(1 - zf 
5(5 — 5z + 5z 2 — z 3 )z 2 \nz 



0.63 



(40) 



15(1 -zf 

where z = m c /mb, and the approximation r(*) ~ z has been made in the order a s corrections. 



The numerical values are for m r 



These values agree with those calculated in Ref. j^j]. 
1.4 GeV and m b = 4.8 GeV. 

The nonperturbative corrections cannot currently be calculated model-independently be- 
cause they depend on the four subleading Isgur-Wise functions that parameterize first-order 
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corrections to the heavy-quark limit, Xi-3( w ) an d v( w )- But they can be estimated using 
QCD sum rules [25] (and, in principle, lattice QCD). In the notation of Neubert 
nonperturbative corrections are 



221, the 



°D* 



>D 



A 



(l/m) 
D* 



A 



(l/m) 
D 



■2 X ;(1)(1 + 2) - 3»(1)(1 - 3-) + 4Vs(l)(l " 3z) 



-j(H-*> 



l-2z + hz 2 
3(1-^) 



7/(1) 



-2.1 



-2(1 + z) [ X , 1 (1)-2 X2 (1) + 6 X ' 3 (1)] + 



2(1 



-277(1) 



1 - 2z + 5z 2 
3(1 -z) 



fa +4 



-1.3, 

8(l-6z + ^)x 2 (i; 
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+ -(1 -3z) X2 (l) -4(1-3^(1)- 



1 + 2 

+ 2(i + *M(i, 9(1 _ 2)2 

r/(l) (5 - 282 + I82 2 - 52,2 s + 252 4 ) 



9(1 



+ 



2t/(1)(1 - 2z + 5z 2 ) (1 + z)(25 - 42,2 + 25z 2 



3(1 



\2 



18(1 - zf 



-2.6p 2 - 1.7 



W(i)- 



-0.3 



1 + 2 



: ] ' 2(1 + z) [ X ;'(l) - 4^(1) + 6xS(l)] - 2r/'(l) ' ] 



1 + 2 



(41) 



where the primes denote d/dw. In these corrections p 2 can be taken to be p 2 DM , since the 
results here do not include corrections of order Q. s Kqp,T>l^ric,b- The parts of these expressions 
for A^» m ^ and A^ m ^ proportional to p 2 disagree with those of Ref. (24}]. 1 The numerical 
estimates are based on the approximate values 7/(1) = 0.6, t/(1) = 0, x'i(l) = 0-3, X2W = 
-0.04, xa(l) = 0.03, and x' 3 (l) = 0.02 |25|. The values t/"(1) = xi'(l) = x 3 '(l) = and 
z = 0.29 were also used. Since these values are model dependent with large uncertainties, the 
numerical estimates of the nonperturbative corrections should be interpreted with caution. 
Reliable lattice results would be a welcome check on such large QCD sum rule estimates of 
these corrections. 

Combining the bounds of Eqs. (|3l)|) and (f3~7|) with Eqs. (J3~9*j) gives the physical bounds 



°DM > 4^D« ( 1 + 



s « + ^ln- m i- 



a 

A 



45tt 4A 2 



13a s , m 2 
ln- 



45tt 4A 2 



5a « A K) , A («.) 



5 A r(l/m) , 

42m c DM + 



* (l/m) 

2mr DM ' 



(42) 



3 / 2 \2 4 2 / 4a s , m 2 



a s („,) 3 A 



(l/m) 



7T 



DW 2 2m, dW 



4a. s m 2 
m- 



45tt 4A 2 



1 The authors of Ref. ,24j have confirmed these findings. They report that the numerical result in their 
Eq. (19) for the difference {a 2 D - cr 2 D ,)/2 changes from 0.17 + 0.20p 2 to 0.17 + 0.29/9 2 . 
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4a, 4a, (cts) a, (cts) 4 A {1/m) A (1/m) 



15tt 5tt " w tt " w 5 2m n UK ' 2m. 
Numerically, these inequalities are 

> \pdA 1 -0.11(0.16)p-0.3 np ] - 0.081(0.009),, + 0.1 

3 / o \2 4 



np j 



<r£, > |(p^j +-p^.[l- 0.066(0.101) p + 0.08 np ] -0.14(0.13) p - 0.003 np , 



> 7 pL[1 + 0.041 (-0.014) p + 0.0„ p ] -0.025(0.0028) p + 0.2 



4 



np ; 



> — (p 2 D ) 2 + -p 2 D [l + 0.025(-0.0089) p + 0.3 np ] - 0.039(0.032) p + 0.1„ p , (44) 

where the values a, = 0.3 (in the MS scheme around 2 GeV) and A = 0.4 GeV have been 
used. The perturbative corrections, with subscript p, are for two values of A. The results 
for A = 2 GeV are first, and those for A = 3 GeV are in parentheses. The nonperturbative 
corrections are labeled by a subscript np. 

Equations (jUj) and fi3|) imply absolute bounds when combined with the corrected form 
of the Uraltsev bound, 



n 2 > 6 + A(Xs In ^ + % (Qs) + A r) {1/m) ( IV 



which comes from Eq. (J30|) and the first of Eqs. (J39)) . and the tree-level Voloshin bound [26], 
p 2 ~ 3/4. A lower bound is required for terms proportional to p 2 D ( t ) with positive coeffi- 
cients, and an upper bound is required for those with negative coefficients. The latter are 
corrections, so the upper bound is only needed at tree level. Note that an upper bound is 
required to estimate the greatest impact the corrections could have on the bound. 
Inserting these inequalities into Eq. (}4*2*|) gives 

2 15 14 «* , K , 3a K) a K) A (1/m) 

> 16 + To7" ln iAi + + ~ A ^> + ^ A ^<*> > ^ 46 ) 

where p|, w is replaced by 3/4 in A^ffi. The absolute bound produced in the same way 
from Eq. (|43|) happens to be identical at leading order and weaker only by the addition of 
— 4a,/157r after perturbative and AQ CD /m c fe corrections are included. Using the numerical 
estimates above, the bounds in Eq. (j4Hj) are 



o 2 D , > 0.94 - 0.26(0.34) p - 0.5 np , 

a 2 D > 0.94 + 0.045(-0.027) p - 0.04 np , (47) 
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where the corrections are labeled as described above. 

When considering the bounds in Eqs. (jUj) and ()46|) and their dependence on A, one must 
bear in mind that the logarithms in the perturbative corrections are only small if A, m^, 
and m c are roughly of the same order. That is the accuracy of the results obtained here. For 
instance, Eq. (|39|) is valid for /i on the order of m c ^, while Eqs. (J36|) and (|37|) are valid for /i 
near A. Therefore, taking the A — > oo limit does not make sense in the absolute bounds. To 
understand the behavior of the bounds in this limit, one would need to sum the logarithms 
of m 2 /A 2 . Since these logarithms are not large for the values of A used here, this extra step 
has been omitted. 

The sum rule bounds derived here should be compared with the dispersive constraints 
usually used to guide the extrapolation of measured form factors to zero recoil. These con- 
straints are derived by computing the vacuum expectation value of a time-ordered product 




_, , , -n 

-0.5 / 


' / 0.5 


1 


1.5 



















FIG. 4: Dispersive constraints on Tjj derivatives combined with the corrected sum rule bounds 
derived here at A = 2 GeV. The interior of the ellipse is the region allowed by the dispersion 
relations. Including the curvature bounds, given by the area above the dashed curves, further 
restricts the allowed region to the shaded area. The darker region is obtained by also including the 
Bjorken and Voloshin bounds. Both perturbative and nonperturbative corrections are included. 
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of b — > c currents in the perturbative regime and then using analyticity to learn about 
the semileptonic regime. The result is equated with a spectral function sum of resonances 
(i.e., a sum of positive quantities). Much as in the derivation of the sum rules here, fo- 
cusing on specific resonances yields form factor constraints. A typical example is shown 
in Fig. HJ The slope and curvature must lie within the ellipse, a constraint that is well 
approximated by a linear relation between the slope and curvature. Data are fit as a func- 
tion of w for \V c b\J r D (*)(l) and p 2 DM , and the second derivative at zero recoil is related to 
the slope according to this relation. For the process B° — > D + £~u, Belle used the relation 
a 2 D /2 = 1.05pi> - 0.15 [9] to find o 2 D = 2.06 ± 0.46 ± 0.29, where the first uncertainty is 
statistical and the second systematic 0]. This value is consistent with the bound above. 

Rather than ^^(w), one typically fits the shape of the axial vector form factor h Al {w), 
which is defined, for example, in Ref. j^J. Like J 7 d*{w), it is equal to the Isgur-Wise function 
in the heavy-quark limit. Its curvature is defined as in Eq. ()38|) and satisfies the bound 
in Eq. ()46|) . with perturbative and nonperturbative corrections given by 



.(«,) 2(l-z)(17-4z+17z 2 ) + 6(9-3z + 4z 2 )zlnz 

0\ = : = 0.65 

Al 27(1 - zf 



S (l/m) = _ 2(1 + z)x[{1) + Azx2{l) + ^ _ + zrj{l) _ +Z ^ 



1+Z 

K) 4(1 - z)(27 - 203z - 68z 2 - 203z 3 + 27z 4 ) - 120(10 + 5z 2 - z 3 )z 2 \nz _ 

a: ~ 225(1 -zf " ° ' 

A2f m) = p 2 [2 zr] (l) -l-z) + 2x?(l)(l +z)- 8z X ' 2 (l) - 4xS(l)(l - 3z) 

+ zt]{1) - 2zr)'(l) - w -0.9p 2 - 0.5 , (48) 

where the numerical estimates are based on the values used above. These values produce 
the absolute bound 

a 2 Al > .94 - 0.071(0.14) p - 0.2 np (49) 

on the curvature of Ha^w). This should be compared with the value found by Belle by 
the procedure described above for the process B° — > D* + e~v. Using the relation <j\ J2 = 
1.08p 2 Ai - 0.23 Q, Belle found a 2 Al = 2.44 ± 0.37 ± 0.41, where the first uncertainty is 
statistical and the second systematic [2^J. This value is also consistent with the bound 
produced here. A plot comparing dispersive constraints and sum rule bounds for this form 
factor would be similar to Fig. but with the sum rule bounds comparatively somewhat 
weaker. 
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V. CONCLUSIONS 



Th is p aper has presented order a s corrections to two new sum rules derived in 



Refs. 



lis pape 

00 



12j in the context of HQET. Section |H] repeated the tree-level derivation 



of a generic sum rule depending on three velocity transfer variables and included one-loop 
corrections. Section ITTT1 studied the axial vector and vector sum rules that result from choos- 
ing specific currents in the generic equation. These led to a s -corrected versions of the sum 
rules of Le Yaouanc et al. for the curvature of the Isgur-Wise function. There were no 
corrections suppressed by the heavy quark masses because the infinite-mass limit was used. 
Section HVI translated these HQET bounds into bounds on physical form factors by including 
perturbative and finite-mass corrections associated with matching HQET onto the full the- 
ory. Numerical estimates were given and compared with experimental values and dispersive 
constraints. The bounds produced here are less powerful than dispersive constraints but 
may provide an important check on those constraints. 



Acknowledgments 

I would like to thank Mark Wise for many discussions and a critical reading of the draft, 
and Zoltan Ligeti for several helpful comments on the draft. This work was supported in 
part by a National Science Foundation Graduate Research Fellowship and the Department 
of Energy under Grant No. DE-FG03-92ER40701. 



[1] E. Eichten and B. Hill, Phys. Lett. B 234, 511 (1990); H. Georgi, Phys. Lett. B 240, 447 
(1990). 

[2] N. Isgur and M. B. Wise, Phys. Lett. B 232, 113 (1989); Phys. Lett. B 237, 527 (1990). 
[3] S. Nussinov and W. Wetzel, Phys. Rev. D 36, 130 (1987). 

[4] M. A. Shifman and M. B. Voloshin, Sov. J. Nucl. Phys. 47, 511 (1988) [Yad. Fiz. 47, 801 
(1988)]. 

[5] M. E. Luke, Phys. Lett. B 252, 447 (1990). 

[6] A. Czarnecki, Phys. Rev. Lett. 76, 4124 (1996); A. Czarnecki and K. Melnikov, Nucl. Phys. 
B 505, 65 (1997). 



18 



[7] S. Stone, in B Decays, 2nd Edition, edited by S. Stone (World Scientific, Singapore 1994), 
p. 283. 

[8] C. G. Boyd, B. Grinstein and R. F. Lebed, Phys. Lett. B 353, 306 (1995); Nucl. Phys. B 461, 
493 (1996); Phys. Rev. D 56, 6895 (1997). 

[9] I. Caprini and M. Neubert, Phys. Lett. B 380, 376 (1996); I. Caprini, L. Lellouch and M. Neu- 
bert, Nucl. Phys. B 530, 153 (1998). 
[10] A. Le Yaouanc, L. Oliver and J. C. Raynal, Phys. Rev. D 67, 114009 (2003). 
[11] A. Le Yaouanc, L. Oliver and J. C. Raynal, Phys. Lett. B 557, 207 (2003). 
[12] A. Le Yaouanc, L. Oliver and J. C. Raynal, Phys. Rev. D 69, 094022 (2004). 
[13] I. Bigi, M. Shifman, N. G. Uraltsev and A. Vainshtein, Phys. Rev. D 52, 196 (1995). 
[14] A. Kapustin, Z. Ligeti, M. B. Wise and B. Grinstein, Phys. Lett. B 375, 327 (1996). 
[15] C. G. Boyd and I. Z. Rothstein, Phys. Lett. B 395, 96 (1997). 

[16] C. G. Boyd, Z. Ligeti, I. Z. Rothstein and M. B. Wise, Phys. Rev. D 55, 3027 (1997). 
[17] N. Uraltsev, Phys. Lett. B 501, 86 (2001). 

[18] E. C. Poggio, H. R. Quinn and S. Weinberg, Phys. Rev. D 13, 1958 (1976). 

[19] A. F. Falk, Nucl. Phys. B 378, 79 (1992). 

[20] N. Isgur and M. B. Wise, Phys. Rev. D 43, 819 (1991). 

[21] J. D. Bjorken, SLAC-PUB-5278 Invited talk given at Les Rencontres de la Valle d'Aoste, La 
Thuile, Italy, Mar 18-24, 1990; J. D. Bjorken, I. Dunietz and J. Taron, Nucl. Phys. B 371, 
111 (1992). 

[22] M. Neubert, Phys. Rept. 245, 259 (1994). 

[23] C. G. Boyd, B. Grinstein and A. V. Manohar, Phys. Rev. D 54, 2081 (1996). 
[24] B. Grinstein and Z. Ligeti, Phys. Lett. B 526, 345 (2002). 

[25] M. Neubert, Phys. Rev. D 46, 3914 (1992); M. Neubert, Z. Ligeti and Y. Nir, Phys. Lett. B 
301, 101 (1993); Phys. Rev. D 47, 5060 (1993); Z. Ligeti, Y. Nir and M. Neubert, Phys. Rev. 
D 49, 1302 (1994). 

[26] M. B. Voloshin, Phys. Rev. D 46, 3062 (1992). 

[27] K. Abe et al. [Belle Collaboration], Phys. Lett. B 526, 258 (2002). 

[28] K. Abe et al. [Belle Collaboration], Phys. Lett. B 526, 247 (2002). 



19 



